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Abstract 

An elementary analytical fluid flow is composed by a geometric do- 
main, a list of analytical constraints and by the function which depends 
on the physical properties, as Reynolds number, of the considered fluid. 
For this object, notions of laminar or weakly turbulent behavior are de- 
scribed using a simple mathematical model. 

1 Elementary analytical fluid flows 

Let L a length unit and a closed domain of M. n , n > 2, formed by a small 
number (e.g. 2 or 3) of L 2 squares or L 3 cubes. Here we consider the plane case, 
but the notions are easly extensible to greater dimensions. Let F a fluid entering 
in Q from a single edge of a square and flowing out from another single edge. Let 
s, < s < 1, a parameter which describes the position of a single fluid particle 
along the inflow edge, so that sqL is the initial position of the generic particle 
associated to the particular value sq of the parameter. Let now s i— » R(s) a 
function R : [0, 1] — > [0, 1] which maps a value of the parameter s to the value 
identifying the position reached by the particle on the outflow edge, so that 
this position is identified by the value R(s)L. If x,y is a cartesian coordinates 
system in the plane of tt, assume that each streamline, or each particle path in 
lagrangian view, is described by a parametric curve t \— > — {x{t), y(t)), with 
< t < 1. The parameter t could not be, in general, the time variable of the 
flow. Let this parametric representation be determined by a set A of analytical 
conditions regarding $(£) and $(£), that is the passage of the particle in some 
suitable points of fl and the velocity field direction in some other (or the same) 
points. 

An elementary analytical flow is a particular list F = {£7, R, A}. We propose a 
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simple mathematical model, at least for a particular example of domain f2, and 
show that the function R(s) can identify the physical properties of the particular 
fluid F flowing in the domain. For simplicity in this work assume that the flow 
is steady. 



2 A mathematical model 

In this section let f2 the domain formed by three L 2 squares, the first two along 
the a;-axis of a cartesian coordinates system from to 2L, and the third above 
the second, from L to 2L y-coordinates. 
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Figure 1: Domain with three unit squares. 



The inflow edge is the segment [(0, 0), (0, L)] of the first square and the outflow 
edge is the segment [(L, 2L), (2L, 2L)\ of the third square. 

Now the set A af analytical conditions is so defined: 

P\. at the inflow edge, for t = and position parameter s, let $(0) = (0, sL); 

Pi. at the outflow edge, for t = 1 and position parameter R = R(s), let 
$(l) = (L(l + i2),2L); 

P 3 . for t — \ a particle path intersects the diagonal line of the second square, 
that is the line of cartesian equation y = —x + 2L\ if a parameter p, < p < 1, 
describes the positions along this line, the condition is 3>(|) = (L(2 — p),pL); 

D]_. at the inflow edge the velocity is parallel to x-axis, so that for every s 
2/(0) = 0; 

D 2 . at the outflow edge the velocity is parallel to y-axis, so that for every s 
x(l) = 0. 

Note that condition P3 is a logical connection between x and y coordinates, so 
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that for each component of <&(i) the independent conditions are four: assuming, 
for semplicity, that the components are polinomial expressions on variable t, 
the candidates to satisfy the set A — {Pi, P2, P3, D\, D2} of conditions are the 
cubics x(t) = at 3 + bt 2 + ct + d and y(t) = et 3 + ft 2 + gt + h. 

Let proceed with the computation of these cubics. 

From derivation condition D\ follows 5 = 0. From passage condition P\ follows 
d = and h = sL. From condition D2 follows c = —3a — 2b. With these partial 
results, from condition P3 follows b = —L(l + R) — 2a and c = a + 2L(1 + R). 
Now, using condition P2, the computation gives a = L(10 — QR — 8p), so that 
the first component is 

x(t) = L(10 - QR - 8p)t 3 + L(-21 + 11R + I6p)t 2 + L(12 - 4R - 8p)t (1) 
With a similar computation, the second component is 

y(t) = L(A + 6s - 8p)t 3 + L(-2 -7s + 8p)t 2 + Ls (2) 
Note that, with this model, the second component doesn't depend on function 

For semplicity, let p — s, that is the particles paths are not disturbed until the 
line y = —x + 2L. In this case the y-component has the form 

y(t) = L{A - 2s)t 3 + L{-2 + s)t 2 + Ls (3) 

3 The laminar case 

In the laminar case assume that particles paths don't intersect themselves, that 
is the outflow positions on the final edge of Cl are the same as inflow positions. 
The analytical expression of this geometrical condition is simply 

R{s) = 1 - s (4) 

The laminar elementary analytical flow F = {17, A, R} is so completely defined; 
its particles paths have the form 

= {L(2 - s){2t 3 - M 2 + 4i), L{2 - s){2t 3 - t 2 ) + Ls} (5) 

Note that the expression Q for R{s) depends on the shape of the domain fj, 
as expected; e.g., in the case of a rectangular domain with the outflow edge 
parallel to the inflow edge, the expression should be R(s) = s. 

4 The weakly turbulent case 

The weak turbulence notion (0) regards the physical situation of a flow where 
turbulence is not fully developed and it could be described by deterministic 



3 



Figure 2: Particles paths for laminar elementary flow. 



mathematical law (@j). For an elementary analytical flow, such mathematical 
description can be made by a suitable choice of the function R(s), which maps 
the initial inflow position of a particle to its outflow position. 
In order to find a possible expression for R(s), assume that some geometrical 
point of the final edge can be reached by more than one fluid particle, so that 
paths intersects themselves, as expected in a turbulent flow. Also, the analyt- 
ical expression of R(s) should be such that, in the case of suitable conditions, 
the laminar case (@J can be a mathematical limit for great values of viscosity 
(viscous limit, see pQ). A mathematical condition is that R(s) £ [0, 1] for every 
s. Assume also that weak turbulence, as expected, is associate to the analytical 
condition R G C m ([0, 1]) for some integer m > 0. 

Note that a possible expression satisfying all the previous conditions is R(s) = 
I [1 + sin(a)], where a — a(s) could be computed using the viscous limit. As- 
sume that this limit is defined for Reynolds number Re approaching a value Cf , 
Re — > cf, with cf a constant which should depend on the nature of the fluid and 
on the geometrical and physical properties of the flow. The most simple form 
for a is a = a Re s + 6, therefore 

Re — > Cf =>■ a — > aefs + b (6) 

Assuming acps + b sufficiently small, for Re w cf follows sin(a) » acys + b, 
therefore it must be 

R(s) = ~ [1 + ac w s + b] = 1 - s (7) 

from which a — — — and 6=1. So the function describing a weak turbulent 
elementary flow is 
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Figure 3: Viscous limit (red) R(s) = | [1 + sin (1 — 2s)], compared to laminar case 
R(s) = l-s. 




Figure 4: Weakly turbulent elementary flows: at left, particle paths for = 3; a t 
right, some particle paths in the case S£ — 2 10 4 . 



5 Concluding remarks and further developments 

The notion of elementary flow is characterized by a domain, a list of analytical 
constraints and by the function R = R(s) : [0, 1] — > [0, 1] which depends on the 
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geometrical and physical properties of the considered fluid. A particular weakly 
turbulent behavior can be obtained when R(s) is a periodic function. Then 
laminar behavior is a viscous limit of this model. 
Some questions arise. 

It could be investigated the relation between particles paths (or streamlines) 
$(£) and the solutions of Navier-Stokes equations, probably obtaining a mathe- 
matical relation between R and pressure, and an expression for the flow constant 

Cg. 

In the case of evolution, e.g. if the t parameter is proportional to time, func- 
tion R(s) could be characterized using the Richardson's Law on particles 
dispersion in turbulent flows. 

Can be full developed turbulence a limit of this weak turbulence model? For 
very high Reynolds numbers the function R(s) expressed by (JSJ formula has no 
limit for every s, or has not a deterministic behavior: the vertical line [0, 1] is 
an accumulation line for its graph, therefore no prediction about the outflow 
position of a particle is possible. Perhaps it can be useful a convergence, for 
some kind of metric, of a succession {R n (s)} n =i,2,3,... of continuous Rs to some 
step function, in order to describe full turbulence flow and its phenomemon of 
intcrmittency ([3]). 

Finally, can be a real fluid flow described by some kind of combination of el- 
ementary flows? Perhaps a combination, or a succession, of elementary flows 
with a common function R(s) and variable analytical constraints {^4 n }n=i,2,3,... 
could be a useful representation of a real flow of a single fluid. 
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